Measurement of the Absolute Branching 
Ratio for the K + — > /i + z/(7) Decay with the 
KLOE Detector. 



The KLOE Collaboration 

O ' 

O ; F. Ambrosino f , A. Antonelli b , M. Antonelli b , C. Bacci^, 

! P. Beltrame b , G. Bencivenni b , S. Bertolucci b , C. Bini-i, 

q ! C. Bloise b , V. Boccij, F. Bossi b , D. Bowring b '°, P. Branching, 

R. Caloi-i, P. Campana b , G. Capon b , T. Capussela f , 
F. Ceradini^, S. Chi b , G. Chiefari f , P. Ciambrone b , 
S. Conetti , E. De Lucia b '\ P. De Simone b , G. De ZorziJ, 
in ■ S. Dell'Agnello b , A. Denig c , A. Di Domenico-i, C. Di Donato f , 

S. Di Falco h , B. Di Mice*/, A. Doria f , M. Dreucci b , G. Felici b , 
A. Ferrari , M. L. Ferrer b , G. Finocchiaro b , C. Forti b , 
P. Franzini j , C. Gatti b , P. Gauzzi^, S. Giovannella b , 
! E. Gorini d , E. Graziani^, M. Incagli h , W. Kluge c , V. Kulikov , 

<d : 
a: 

; M. Moulson b , S. Miiller c , F. Murtas b , M. Napolitano f , 

F. Nguyen^, M. Palutan b , E. Pasqualucci •' , A. Passed e , 



(N 



X 



F. Lacava^, G. Lanfranchi b , J. Lee-Franzini b ' m , D. Leone c , 
M. Martini b , P. Massarotti f , W. Mei b , S. Meola f , S. Miscetti b , 



c 



V. Patera 13 ' 1 , F. Perfetto f , L. Pontecorvo-i, M. Primavera d , 
P. Santangelo b , E. Santovetti k , G. Saracino f , B. Sciascia b , 
A. Sciubba b '\ F. Scuri h , I. Sfiligoi b , T. Spadaro b , M. Testa 
L. Tortora^, P. Valente-', B. Valerianic G. Venanzoni b , 
S. Veneziano^, A. Ventura d , R.Versaci ' 1 , G. Xu b ' a 

^Permanent address: Institute of High Energy Physics of Academica Sinica, 

Beijing, China. 

h Laboratori Nazionali di Frascati delVINFN, Frascati, Italy. 

c Institut fur Experimentelle Kernphysik, Universitdt Karlsruhe, Germany. 

d Dipartimento di Fisica delVUniversita e Sezione INFN, Lecce, Italy. 

Permanent address: Institute for Theoretical and Experimental Physics, Moscow, 

Russia. 

f Dipartimento di Scienze Fisiche delVUniversita "Federico II" e Sezione INFN, 



Preprint submitted to Elsevier Science 



7 February 2008 



Napoli, Italy 

^•Permanent address: Budker Institute of Nuclear Physics, Novosibirsk, Russia. 

h Dipartimento di Fisica dell'Universita e Sezione INFN, Pisa, Italy. 

1 Dipartimento di Energetica dell'Universita "La Sapienza" , Roma, Italy. 

J Dipartimento di Fisica dell'Universita "La Sapienza" e Sezione INFN, Roma, 

Italy. 

k Dipartimento di Fisica dell'Universita "Tor Vergata" e Sezione INFN, Roma, 

Italy. 

Dipartimento di Fisica dell'Universita "Roma Tre" e Sezione INFN, Roma, Italy. 

m Physics Department, State University of New York at Stony Brook, USA. 

n Permanent address: High Energy Physics Institute, Tbilisi State University, 

Tbilisi, Georgia. 

° Physics Department, University of Virginia, USA. 



Abstract 

We have measured the fully inclusive K + — > fi + u{^) absolute branching ratio 
with the KLOE experiment at DA$NE, the Frascati ^-factory. From some 865,283 
K + — > /x + i/(7) decays obtained from a sample of ~5.2 x 10 8 0-meson decays, we 
find BR(A + — > ^ + f M (7))=0.6366±0.0009 sta t. ±0.0015 syst ., corresponding to an over- 
all fractional error of 0.27%. Using recent lattice results on the decay constants of 
pseudoscalar mesons one can obtain an estimate for the CKM mixing matrix ele- 
ment |KJ=0.2223±0.0026. 



1 Introduction 

The most recent measurement of the K — > \xv branching ratio, ref. 1, based 
on 62,000 events, dates back to 1972, more than 30 years ago, and relies on a 
sample of ~10 5 kaon decays. The authors of ref. 1 quote an error of ~0.7%, 
the statistical error due to the event count being 0.4%. This error is reduced 
in the PDG fit [2] to 0.27% and the value changed by 0.3%. While the proce- 
dure is correct in principle, it can lead to incorrect results because of incorrect 
data included in the fit. This has been the case sometimes, leading among 
other things to the suggestion, also in ref. 2, of a value for |T4s| apparently 
inconsistent with unitarity of the CKM mixing matrix [3] . Another problem 
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with measurements performed more than 30 years ago is due to the fact that 
the effect of radiative corrections was not fully appreciated. It is therefore 
impossible to understand what fraction of the radiative decay is included in 
the quoted results. Inclusion of all radiation is however necessary to compare 
to models or to extract fundamental parameters such as a coupling constant. 
With all the above in mind as well as recent developments in numerical or 
lattice QCD calculation, we have begun a program of new precise, fully inclu- 
sive, kaon branching ratio measurement. 

We report in the following a measurement of BK(K + — > fi + u('~f)) performed 
with the KLOE detector at DA$NE, the Frascati 0-factory. The measure- 
ment is based on an integrated luminosity of ~175 pb _1 , collected in 2001-02. 
DA$NE is an e + e~ collider operated at a total energy of W — 1020 MeV, the 
mass of the 0(lO2O)-meson. Equal energy positron and electron beams collide 
at an angle of (n— )25 mrad and produce 0-mesons with a small transverse 
momentum of ~12.5 MeV/c. The collision frame moves therefore in the lab- 
oratory with a velocity /3~0.0125. In the center of mass, the 0-meson decays 
into anti-collinear K pairs of ~125 MeV/c momentum. In the laboratory 
this remains approximately true: detection of a K ± tags the presence of a 
K T of given momentum and direction. The decay products of the K ± pair 
define two spatially well separated regions called in the following the tag and 
the signal hemispheres. Identified K T decays tag a beam and provide an 
absolute count. This procedure is a unique feature of a ^-factory and provides 
the means for measurements of absolute branching ratios, i.e. ratios Tj/T^ 
rather than ratios of BR's Ti/Tj. 



2 The KLOE detector 

The KLOE detector consists of a large volume drift chamber and a sampling 
calorimeter. The drift chamber (DC) [4], of 3.3 m length and 2 m radius, has 
a full stereo geometry and operates with a 90% helium-10% isobutane gas 
mixture. Tracking in the DC provides measurements of the vector momentum 
of charged particles with cr(p±)/pi_ < 0.4% and two track vertices to 3 mm. In 
the following we use a coordinate system with the z-axis defined as the bisec- 
trix of the e + e~ beams, the y-axis vertical and the x-axis toward the center 
of the collider rings. 

The calorimeter (EMC) [5] consists of a cylindrical barrel and two endcaps cov- 
ering a solid angle of 98% of Aw. Photons showering in the lead-scintillator-fiber 
EMC structure are detected as local energy deposits by clustering signals from 
read-out elements. The calorimeter information consists of energy, position of 



entry point and time of arrival with accuracies of (Je/E = 5.7%/ J E(GeV), 
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a z = 1.2cm/ ^(GeV), = 1.2 cm and a t = 54 ps/y^GeV) © 50 ps. 
Calorimeter clusters not associated with a DC track indicate arrival of neu- 
tral particles and the computed time of flight identifies photons with excellent 
precision. Time of flight also allows good separation of electrons from muons, 
pions and kaons. 

A superconducting coil surrounds the entire detector and produces a solenoidal 
field B = 0.52 T. The trigger [6] is based on the detection of isolated energy 
deposits in the calorimeter and on hit multiplicity in the drift chamber. Only 
events triggered by the calorimeter have been used in the present analysis. 
This choice ensures a far more reliable estimate of all necessary efficiencies. 
The trigger system also includes a veto for cosmic-ray muons (cosmic ray veto 
or CRV) based on energy deposits in the outermost layers of the calorimeter 
and followed by a third-level software trigger able to identify most of the 
events. A software filter (filfo or FLF), based on the topology and multiplicity 
of calorimeter clusters and drift chamber hits, is applied to filter out machine 
background. 

Both CRV and FLF are sources of event rejection. Their effect on the BR 
measurement has been studied on control samples which do not undergo, re- 
spectively, the CRV, and the FLF filter. 



3 The measurement 

The entire data sample of 175 pb -1 is divided into two subsamples. 
Some 60 pb -1 of data have been used for the BR measurement. The remain- 
ing 115 pb -1 have been used to evaluate efficiencies and the background. The 
branching ratio measurement [7] is based on the use of K~ — > \x~v decays for 
event tagging and to search for the K + — > fi + u( , -f) signal among K + decays. 
The tagging selection is based on the presence of a two-tracks vertex in the 
DC which signals the K~ decay [7]. 

Nuclear interactions, NI, of kaons affects the branching ratio measurement 
but not the tagging procedure. Since <T N i(K + )~a m (K~)/10 2 , the choice above 
minimizes the corrections to be calculated to account for the effect. The cor- 
rections are in fact negligible. The large number of K + decays is sufficient for 
achieving a statistic accuracy at the 0.1% level, comparable with the system- 
atic error. 

To avoid any bias due to differences in the trigger efficiency among the K + 
decay modes on the signal "hemisphere" , the particles on the tagging side are 
required to deposit enough energy in the calorimeter to trigger the data acqui- 
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sition. Nevertheless, the tagging criteria exhibit a residual, small dependence 
on the decay mode of the K + on the signal hemisphere, introducing a tag bias 
or TB, that has been studied using Monte Carlo simulation (MC) samples [8] 
and checked on data. 

The search for positive kaon moving outwards in the DC, with momentum 
70 < pk < 130 Mev/c, is performed on the sample of tagged events. The 
point of closest approach to the beam line with coordinates {x, y, z} is eval- 
uated extrapolating the kaon track backwards to the beam line, taking into 
account the kaon energy loss. Kaon tracks with \z\ < 20 cm and y/x 2 + y' 2 <10 

cm, and kaon decay vertices in the fiducial volume, 40 < \/x v + y v < 150 cm, 
are selected. 

The number of K + — > /i + z/(7) decays is obtained counting the events with 
225< p* <400 MeV/c; p* is the charged decay particle momentum computed 
in the kaon rest frame assuming the pion mass. The p* distribution is shown 
in fig. 1. 




160 180 200 220 240 260 280 300 



Fig. 1. Monte Carlo spectra of the charged decay particle momentum transformed 
to the kaon rest frame, assuming the pion mass. The two peaks corresponds to pions 
from K + -> tt+tt at 205 MeV/c and K + -> y+u at 236 MeV. The black/white line 
represents the signal while the grey line the signal plus background. Contributions 
from K + — > y + u^f are also shown. 

The spectrum in fig.l, obtained from MC simulation, shows a 2% contami- 



5 



nation from various background sources, namely K + — > 7r + 7r°, K + — > tt 1 + v. 
Since the maximum momentum of the pions from three-pion decays is 125 
MeV/c, these channels do not contribute to the background. All the back- 
ground sources in this analysis have one neutral pion in the final state. The 
neutral pions are identified by detecting the photons from tt° — > 77 decay. 
The photons are identified as isolated energy deposits in EMC not associated 
with tracks and satisfying the constraints on 7r° mass reconstruction and time 
correlation with the kaon decay vertex [7]. This selection allows us to obtain 
directly from data the p* distribution of the background. 

The p* distribution for the signal events has been obtained from a data control 
sample described in the following efficiency evaluation. This distribution has 
been used together with the shape of the background sources to fit the overall 
p* spectrum, figure 2 left, and to perform background subtraction. The result, 
after background subtraction, is shown in figure 2 right. 




) r ... 1 ... 1 ... 1 ... 1 ... 1 ... ii I 

180 200 220 240 260 280 300 220 260 300 340 380 



Fig. 2. Left: Spectrum of the charged-secondary momentum in the kaon rest frame 
for data. The solid line indicates the fit done with a linear combination of signal 
and background distributions obtained from data. Right: The same spectrum af- 
ter background subtraction. The range used for counting K + — > /x + i/(7) events is 
indicated. 



The branching ratio is obtained from: 

BR(K + -> „+„(7)) = N *^M) x 1 (1) 

iv Tag e ^CRV ^FF ^TB 



where N^+^ the events selected as signal (hereafter signal count), 

N Tag is the number of tagged events and e is the efficiency Cqrv and Cff are 
the corrections for the effects due to the cosmic-ray veto and the filfo proce- 
dure, and Ctb accounts for the tag bias effect. 
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The efficiency of the analysis cuts has been determined directly on data using 
a control sample of K — > ^(7) events selected exploiting their typical signa- 
ture in the EMC. The control sample consists of events with K~ — > /i~i/(7), 
providing the tag, and signal events K + — > /i + z/(7) selected using only EMC 
information. This criterion is mostly independent from the selection procedure 
based on DC information that has been used for obtaining the signal count. 
The EMC selection requires only one cluster with E > 80 MeV (High Energy 
Cut), plus any number of clusters with energy below 20 MeV (Low Energy 
Cut) which can be due to photons from K + — > fi + w-f . Further cuts on the 
energy (E c {) and on the distance in the tranverse plane from the z-axis (R xy ) 
of the cluster are applied to get rid of machine background and spurious clus- 
ters. Namely we require the cluster either to be on the barrel R xy >197 cm or 
to satisfy the relation E c \/ '(lGeV) + R xy /b(cm) > 110. A correction of about 
0.1%, due to a tiny difference between the efficiency evaluated on the control 
sample and the selection efficiency on the signal sample, has been estimated 
from MC. The efficiency is e = 0.3153 ± 0.0002. 

The corrections to the event rejection described above, Ccrv=1-0005 and 
Cpp — 1 = O(10~ 5 ), have been directly measured on control samples which 
do not undergo, respectively, the cosmic-ray veto, and the filfo filter. The cor- 
rection for the tag bias, Ctb=1-0164±0.0002, has been evaluated on Monte 
Carlo samples and the distribution of the variables used for the tag selection 
have been checked on data. The following sources of systematic uncertainties 
have been studied varying the selection cuts: 

- the requirements on the tagging hemisphere; 

- the trigger requirements; 

- the definition of the fiducial volume; 

- the background evaluation procedure; 

- the choice of the [p* minJ p* max ] range; 

- the energy cuts for the efficiency sample; 

- the effect of high-energy radiative photons (E 1 >20 MeV). 

Furthermore, the dependence of the measurement on the charged kaon life- 
time, which affects the estimate of the geometrical acceptance, has been stud- 
ied varying the lifetime value used in the MC simulation. The effects due to K + 
nuclear interactions have been evaluated from MC simulation and measure- 
ments available in literature. The stability of the measurement with respect 
to different data taking conditions has been checked. The corresponding sys- 
tematic errors are listed in table 1. 

The statistical error due to the event count is 6 x 10~ 4 and becomes 9 x 10~ 4 
including the statistics of MC simulation and data used for the efficiency eval- 
uation. 
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oource 


Value 


Low Energy Cut 


5 x 1U 


bj^ >2() MeV 


7 x 10 


High Energy Cut 


o v 1 n — 4 
z x 11) 


Fiducial Volume 


5 x ID 


Background 


o \ / in — 4 

o x 10 


p* range 


O X 11) 


lag definition 


1 X 11) 


ML Lifetime 


< ID 


Nuclear interactions 


A vy 1 n — 4 

< 4 X 1U 


r llio 


< o X 1U 


Cosmic ray veto 


O(10" 6 ) 


Trigger 


9 x 10" 4 


Total syst. 


15 x icr 4 



Table 1 

Summary table of systematic uncertainties. 
4 Conclusions 

BR Measurement 

On a sample of tagged events Nrags = 4, 237, 329, we found a number of signal 
events, with 225 < p* < 400 MeV/c, of N K +^+ u{l) = 865, 283. 
Using eq. 1, the absolute branching ratio is: 

BR(K + /i + z/(7)) = 0.6366 ± 0.0009 stat ± 0.0015,^. (2) 

The KLOE measurement of the BK(K + —> /i + z/(7)), fully inclusive of final- 
state radiation, has a 0.27% accuracy and it is based on an unprecedented 
statistics and carefully controlled systematics. 

V us extraction 

The recent publication [9] of the results of lattice QCD calculations has re- 
newed the interest in improving the accuracy of the BK(K + — > /i + z/(7)), which 
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represents an experimental alternative to the semileptonic kaon decays in mea- 
suring \V US \ as pointed out by Marciano in ref. 10. The extraction of this CKM 
matrix element is based on the ratio of the decay rates for the inclusive decays 
K + — > n + v{p() and ix + — > /i + z/(7): 

T(K -> Hj) = mK { l ~ ^) \V US \ 2 fl_ 1 + fCg 

r(7r^M7))" m \Vu d \ 2 ft i + [) 



where fx and /„. are, respectively, the kaon and the pion decay constants; 
C n and Ck parametrize the radiative-inclusive electroweak corrections, taking 
into account bremsstrahlung emission of real photons and of virtual-photon 
loop contributions as well. Using the branching ratios of K + — * /i + i/(7) and 
7r + — ► n + v{^() decays, the \V u d\ value from super-allowed nuclear beta decays, 
C n and Ck from ref. 10 and references therein, and the new lattice calculation 
of fx/ fn [9], it is possible to extract |T4s| with an uncertainty at the percent 
level, whose error is mainly dominated by the accuracy of lattice calculations. 
From the BR measurement, using the determination of fx/fir — 1.210(4)(13), 
we have obtained the ratio: 



V us 



V, 



ad 



2 

= 0.05211 ± 0.00016 ± 0.00019 ± 0.00117 (4) 



where the errors correspond, respectively, to the experimental, the structure- 
dependent radiative corrections, and the lattice uncertainties. Taking V u d from 
super-allowed nuclear (3 decays, V u d = 0.9740 ± 0.0005 [11], one determines 
\V US \: 

\V us \ K +^+u(-y) = 0.2223 ± 0.0026. (5) 



The accuracy is dominated by the knowledge of fx / fn from lattice calculation. 
Alternatively, the unitarity relationship \V u d\ 2 = 1 — | | 2 can be assumed in 
eq. 4 giving: 

| V us Unitarity = 0.2225 ± 0.0025 (6) 



independent of the V u d measurement. The result quoted in eq. 5 is in agreement 
with the unitarity of the CKM matrix (eq. 6) and with the determinations 
of V us from semileptonic decays, whose precision is dominated by the /+(0) 
calculations [12]. 
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